arXiv:astro-ph/0211390vl 18 Nov 2002 


1 


Algebraic Properties of the Real Quintic Equation for a Binary 

Gravitational Lens 

Hideki AsADA^’^Ji Taketoshi Kasai^ and Masumi Kasai^’0 

^ Faculty of Science and Technology, Hirosaki University, Hirosaki 036-8561, Japan 
^ Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1 
D-85741 Garching, Germany 

(Received ) 

It has been recently shown that the lens equation for a binary gravitational lens, which is 
apparently a coupled system, can be reduced to a real fifth-order (quintic) algebraic equation. 
Some algebraic properties of the real quintic equation are revealed. We find that the number 
of images on each side of the separation axis is independent of the mass ratio and separation 
unless the source crosses the caustics. Furthermore, the discriminant of the quintic equation 
enables us to study changes in the number of solutions, namely in the number of images. 

It is shown that this discriminant can be factorized into two parts: One represents the 
condition that the lens equation can be reduced to a single quintic equation, while the other 
corresponds to the caustics. 


§1. Introduction 


Gravitational leasing by binary objects has attracted a great deal of interest 0^ 
since the pioneering work of Schneider and Weifi (1986). cP It offers an important 
tool in the search for extrasolar planets, i).0) which is indepernlent of the Doppler 
method, which suffers from uncertainties in inclination angles. B* ** However, it is not 
easy to study binary leasing in detail: The lens equation adopted until now consists 
of a set of two coupled real fifth-order (quintic) equations, ol muivalently, a complex 
quintic equation0 that is expressed in a complex notation. As a result, it seems 
quite difficult to solve these equations, since there are no well-established method 
for solving coupled nonlinear equations numerically with sufficient accuracy. B 

For a symmetric binary with two equal masses, the lens equation for a source on 
the symmetry axes of the binary becomes so simple that we can find analytic solu¬ 
tions.™ For star-planet systems, furthermore, the ma^ratio of the binary is so small 
that we can find approximate solutions in general. Such approximate solutions 

have been us^ to study the shift of the photocenter position due to astrometric 
microlensing. EdP 

One of the present authors recently showed that the lens equation can be reduced 

S a single master equation that is fifth-order in a real variable with real coefficients. 

As a consequence, the new formalism provides a much simpler method to study 
binary leasing. In this paper, we reveal algebraic properties of the real quintic 
equation that have not yet been studied.The discriminant of the quintic equation 
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is shown to be factorizable into two parts: One represents the condition that the lens 
equation can be reduced to a single quintic equation, while the other corresponds to 
the caustics. 

The caustics are curves on which the Jacobian of the lens mapping vanishes on a 
source plane. The number of images has been classified in terms of caustics. Caustics 
for two-point masses have been investigated in detail on the basis of a set of two 
coupled real fifth-order equations as an application of catastrophe theoryEu^ and with 
a complex formalism. US In binary leasing, three images appear for a source outside 
the caustics, while five images appear for a source inside the caustics. Regarding the 
number of images, it should be straightforward to count the number of solutions for 
the lens equation. However, this has not been done for a binary lens system, because 
the lens equations used are coupled nonlinearly. One of the main purposes of the 
present paper is to count these solutions. 

This paper is organized as follows. First, the derivation of the real quintic lens 
equation is briefly summarized. Then, a condition for its reduction is presented. 
Next, the discriminant of the equation is computed, and in this way we clarify the 
relation between the discriminant and the caustics. 

§2. Discriminant of a real quintic lens equation 
2.1. A real quintic lens equation 

We consider a binary system of two point masses Mi and M 2 and a separation 
vector L from object 1 to 2, which is located a distance from the observer. In 
units of the Einstein ring radius angle 0-^, the lens equation reads 

= + ( 2 - 1 ) 


where j3 and 6 denote the vectors representing the positions of the source and image, 
respectively, and we defined the mass ratios ui and V 2 and the angular separation 
vector I as 


Ml 

“ Ml + M 2 ’ 

(2-2) 

M 2 

~ Ml + M 2 ’ 

(2-3) 

<= ^ , 

(2-4) 

We have the identity z^i -l-z .^2 = 1- For brevity, U 2 is denoted by n. 

Equation (|2T|) is a 


set of two coupled real quintic equations for {0x,9y), equivalent to a single complex 
quintic equation for 9^ + 

Let us introduce polar coordinates whose origin is located at the position of ob¬ 
ject 1 and whose angle is measured from the separation axis of the binary. The coor¬ 
dinates for the source, image and separation vector are {Px,ldy) = (pcos psin y?), 
{9x,9y) = (rcos (/), rsin(/)) and {ix,^y) = (f’, 0), respectively, where p, r and f > 0. 
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For brevity, cos (p and sin (p are written as C and S', respectively. The parallel and 

(2-5) 


vertical parts of Eq. (^) with respect to the separation vector i are 


1 — V 

r cos 0(1 - 2 — 


r sin( 


1 - 


r^ 
l-v 


r^ — 2ir cos 0 + 

V 

— 2ir cos (fpf? 


= pC- 
= pS. 




r 2 _ 2£r cos 0 + f^ ’ 


( 2 - 6 ) 


Two particular cases p = 0 (for which the source is located behind the lens object 
Ml) and S = 0 (for which the source is located on the axis) have been studied, 0 
where analytic solutions for an arbitrary mass ratio were obtained explicitly. In the 
following, we concentrate^ on the case p 7 ^ 0 and S 7 ^ 0 (i.e. off-axis sources). In this 
case, it has been shown 113) that Eq. (^) reduces to a fifth-order equation for cot 0, 


^aj(cot 0 )® * = 0 , 


i=0 


where we have defined 


(2-7) 


(2-8) 


ao = uip^S^, 

ai = ^2^2 ^ 2VC52 _ £2(2^2^2 _ ^ 4 ^ 2 ^ _ 2fp^CS^ + p'^S'^ 

-u{5£p^CS^ - Afp'^S^), (2-9) 

02 = -2p'^CS - Up^iC^S - S^) + f{‘ip‘^CS - 2p^C5) + 4e^p^C^S - 2i^p^CS 


+v 


i{2ps + ^p^c'^s - 2p'^s^) - iQep^cs + 2eps 


.3 c3 


2 „2/ 


( 2 - 10 ) 


03 = p2 + ^(2p3^3 _ 10^3^52) ^ £2(_2^2^2 ^ ^ ^4) _ ^ ^4^; 


i('3 


oi „2 




i{-2pC - + 6p-^CS^) + 6tp^C^ - 2tpC 


+ 12^ f, 


04 = -2p^CS + 8£p^C^S - t{4.p^CS + 2p^CS) + U^p^C'^S - 2tp^CS 


3 „3/^2 I 


fA „ 2 / 


+V 


£{2pS - 4.p'^C^S + p'^S^) - 2tp^CS + 2tpS 


05 = p^C^ - 2£p^C'^ + e{2p^C'^ + pV2) - 2^^p^C^ + tp^C 


.3^3 


„ 2^2 


4^2\ 


/j3 ^3/^3 


«4 „2r<2 


+ 12 


-£{2pC + p-^CS^) + 2tp^C^ - 2tpC 


+ 12^ 


2o2 


( 2 - 11 ) 


( 2 - 12 ) 


(2-13) 


Here, none of these coefficients oq, ■ • • 05 are singular, because they are polynomials 
in £, p, Z 2 , C and S all of which are finite. 

Equation (2-7) for cot0 has an advantage over the original formulation 0 for 
tan 0 , 

5 

^ Oi(tan 0 )* = 0 . ( 2 T 4 ) 

i=0 

In particular, because 05 , the coefficient of tan^ 0, can vanish, tan0 can be infinite, 
while, as seen from Eq. (^), cot 0 remains finite, because oq never vanishes for 
off-axis sources. Hence, any solution of Eq. (2-7) is finite. As a consequence, we 




















4 


H. Asada, T. Kasai and M. Kasai 


can derive the following theorem: When the source moves around off the symmetry 
axis {Py = 0, i.e., S = 0), images never cross the separation axis {By = 0, i.e., 
cot (j) = ±oo). This implies that the number of images on each side of the separation 
axis is invariant, unless the source crosses the caustics. This property holds for 
arbitrary I and v, because themumber of images is an integer. An analytic study of 
a symmetric binary {v = 1/2)B1’113 is sufficient to determine the invariant number. 
In this case, only a single image appears on the source side of the axis, while two 
images appear on the side opposite to the source outside the caustics or four images 
appear on the side opposite to the source inside the caustics. 

We also hnd the relation involving 


where we defined 


r sin (j) = 


lTo(cot (/>) 
Wi (cot p) ’ 


(2-15) 


Wo (cot p) = {tpC - ip^C^ -ni + pC) 

—pS{f‘ — 2ipC + 1) cot p — ip'^S'^ cot^ 4>, 


lTi(cot (j)) = p{C — S cot (j)) 
This can be recast into the form 


pS cot^ (j) + 2{l — pC) cot 4> — pS 


(2-16) 

(2-17) 


r = n{(l),p,ip,£,u), 


where 77 is a function of p, p, (p, £ and n. 

It should be noted that Eq. (2T5) is obtained from the two equations 


— 2£r cos p + = 


ui 


p{C — S cot (p) ’ 


i — p{C — S cot (p) 


sin (p — ipSr — £{1 — u) sin </> = 0, 


(2-18) 


(2-19) 

( 2 - 20 ) 


where ITi 7 ^ 0 has been assumed. Let us consider the case Wi = 0 in detail. This 
means that Wq = 0 also. Substituting Eq. (2-19) into in Eq. (2-20) gives us 
the equation Wq = 0 for cot (p, because the coefficients of r vanish, due to Wi = 0 . 
For Wo = Wi = 0, Eqs. (2T9) and (2-20) become identical. Now, we determine 
the condition for which there exists a^t cp satisfying the two equations Wq = 0 and 
Wi = 0. It is given by the resultantli^P defined for Wq = uq cot^ (p + ui cot cp + U 2 
and Wi = Vo cot^ cp + vi cot^ cp + V 2 cot (p + as the determinant 


Uo 

Ul 

U2 

0 

0 

0 

Uo 

Ul 

U2 

0 

0 

0 

Uo 

Ul 

U2 


Vi 

V2 

V3 

0 

0 

Vo 

Vi 

V2 

V 3 


Using Mathematica, this is evaluated as 


( 2 - 21 ) 


W = S^p'^ 


72/ _ 4Cfp^ - p2(l _ £4 _ 4(^2^4 _ 
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+2C^p{l - - fv{2 - 21^ - v) 


( 2 - 22 ) 


If and only if II^ = 0, there exists a cot (j) satisfying Wq = Wi = 0. In this case, 
we must solve Wq = Wi = 0 for cot f. Then, substitution of the solution into Eq. 
(2T9), which is identical to Eq. (2-20), gives us a quadratic equation in r. This can 
be solved immediately. This yields two solutions on a line dehned hy y/x = coif). 
We should note that VE = 0 is an artifact of our polar coordinates. In the following, 
we consider only the case W 0. 


2.2. Discriminant 

We factorized the discriminant 0 of the quintic equation Eq. (^) using 
Mathematica0 into the form 




(2-23) 


where we defined 


16 


K = ^kip 

i=0 


16-i 


(2-24) 


for ko = ki = SCf{l — — 2i/), • ■ ■. Here, A: 2 , • • •, ki% have rather lengthy forms. 

Changes in the number of images are usually studied in terms of the caustics 
(J = 0). In our approach for a real algebraic equation, changes in the number of 
solutions for the lens equation Eq. o correspond directly to those in the signature 
of the discriminant D^. In the next subsection, we show that = 0 is a condition 
for the caustics. 


2.3. Caustics 

Caustics are curves on the source plane, where the Jacobian of the lens mapping 
\d(3/dd\ vanishes. In our polar coordinates, the Jacobian is expressed as 


J = 


(9/3 


P 1 

r d{'r,4>) ' 
d{p,>p) 


(2-25) 


Since r and p are nonvanishing and finite, \dl3/d6\ vanishes where 


d{r, f) 

d{p,ip) 


(2-26) 


Differentiation of Eq. 


(2-7) gives 


<Pd(l) 


5 


i=0 


ddi , do^i , 1 „ 

= 0. 

op op 


Here, we dehned 


= '^bi{cot f) 


6—z 


2=0 


(2-27) 


(2-28) 
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where 6 o = 5ao, bi = 4ai, 62 = 3 a 2 + 5ao, 63 = 2 a 3 + 4ai, 64 = 04 + 802 , 65 = 203 
and he = 04 - Therefore, we find 


9/3 ^ 9p 


i=0 

5 




The total derivative of Equation (|2T8[) leads to 

/ flK flK,^ 

^dr = 


BU dn^. 


dp 


Bn Bn^. ^^,5_iBai 


d</9, 


where we used Eq. ( p-27| ) to eliminate dcf). Hence, we obtain 

Br 
Bp 
Br 
Bip 

so that we find the Jacobian in the polar coordinates as 

B{r, 4>) Bn B4> Bn B(j) 

B{p, ip) Bp Bip Bip Bp ’ 


Bn Bn B(f) 
Bp B4> Bp^ 
Bn Bn Bcj) 
Bp B(j) Bp' 


(2-29) 

(2-30) 


(2-31) 

(2-32) 

(2-33) 

(2-34) 


where Bn/Bp and Bn/Bp are finite. The Jacobian B{r,4>)/B{p,p) is infinite if and 
only if (^ = 0 , because cotcp, Bai/Bp and BaijBp in B(j)/Bp and B(j)jBp are all finite. 

The lens equation and the Jacobian set equr^o 0 are considered a set of equa¬ 
tions for cot (j)- The vanishing of the resultantllBl is the necessary and sufficient 
condition that there exists a cot (j) satisfying both equations. The resultant becomes 


aO 

al 

a 2 

a3 

a4 

a5 

0 

0 

0 

0 

0 

0 

aO 

al 

a 2 

a3 

a4 

a5 

0 

0 

0 

0 

0 

0 

aO 

al 

a 2 

a3 

a4 

a5 

0 

0 

0 

0 

0 

0 

aO 

al 

a 2 

a3 

a4 

a5 

0 

0 

0 

0 

0 

0 

aO 

al 

a 2 

a3 

a4 

a5 

0 

0 

0 

0 

0 

0 

aO 

al 

a 2 

a3 

a4 

a5 

60 

61 

62 

63 

64 

65 

66 

0 

0 

0 

0 

0 

60 

61 

62 

63 

64 

65 

66 

0 

0 

0 

0 

0 

60 

61 

62 

63 

64 

65 

66 

0 

0 

0 

0 

0 

60 

61 

62 

63 

64 

65 

66 

0 

0 

0 

0 

0 

60 

61 

62 

63 

64 

65 

66 


which is a quite complicated function of u, i, p, C and S. Using Mathematica, tD 
we factorized it into 


J 5 = -6452^^2(1 - - pf + ^ 2 ^^ 


(2-36) 
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where {pC — is always positive for off-axis sources {S 7 ^ 0). Hence, from 

Eqs. (|2^ and (l2^ , we have shown that = 0 corresponds to the caustics 
(J 5 = 0 ). Actually, A' = 0 is equivalent to the equation for the caustics obtained 
using the complex formalism. ’113 


§3. Conclusion 

Some algebraic properties of the real quintic equation for a binary gravitational 
lens have been studied. We find that the number of images on each side of the 
separation axis is independent of the mass ratio and separation unless the source 
crosses the caustics. This result should be useful for checking numerical results. 
Furthermore, the discriminant of the quintic equation was shown to factorize into 
two parts: One represents the condition that the lens equation can be reduced to a 
single quintic equation, while the other corresponds to the caustics. 
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